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Abstract 

Perturbative Quantum Chromodynamics (pQCD) predicts that the smah-x gluons in a hadron 
wavefunction should form a Color Glass Condensate (CGC), characterized by a saturation scale 
Qs{x,A) which is energy and atomic number dependent. In this paper we study the predictions 
of CGC physics for electron - ion collisions at high energies. We consider that the nucleus at high 
energies acts as an amplifier of the physics of high parton densities and estimate the nuclear struc- 
ture function F2'{x, Q"^), as well as the longitudinal and charm contributions, using a generalization 
for nuclear targets of the lancu-Itakura-Munier model which describes the ep HERA data quite 
well. Moreover, we investigate the behavior of the logarithmic slopes of the total and longitudinal 
structure functions in the kinematical region of the future electron - ion collider eRHIC. 
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In the high energy hmit, perturbative Quantum Chromo dynamics (pQCD) predicts that 
the small-x gluons in a hadron wavefunction should form a Color Glass Condensate (CGC), 
which is described by an infinite hierarchy of coupled evolution equations for the correlators 
of Wilson lines ^^SSQl- In the absence of correlations, the first equation in the Balitsky- 
JIMWLK hierarchy decouples and is then equivalent to the equation derived independently 
by Kovchegov within the dipole formalism Q]. The Balitsky-Kovchegov (BK) equation 
describes the energy evolution of the dipole-target scattering amplitude A/'(x, r). Although a 
complete analytical solution is still lacking, its main properties are known (for recent reviews 



see, e.g. 



analytical s( 

mm 



(a) for the interaction of a small dipole (r ^ 1/Qs), A/'(r) 



implying that this system is weakly interacting; (b) for a large dipole (r ^ 1/Qs), the 
system is strongly absorbed and therefore //{r) ^ 1. This property is associated to the 
large density of saturated gluons in the hadron wave function. Furthermore, several groups 
have studied the numerical solution of the BK equation lOL 111, 
the theoretical predictions. In particular, the studies presented in |ll|, 
that the BK solution for fixed constant coupling preserves the atomic number dependence 



and confirmed many of 
, Q have demonstrated 



of the saturation scale present in the initial condition, while for running ag this dependence 



is reduced with increasing rapidity, as predicted by Mueller in Ref. 
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The search for signatures of parton saturation effects has been an active subject of research 



in the last years (for recent reviews see, e.g. 



QflQ). On 



one hand, it has been observed 



that the HERA data at small x and low can be successfully described with the help of 



saturation models [1^ lla, llll, llSl , with the experimental results for the total cross section 
and inclusive charm production |2l| presentinsr the property of geometric scaling. On 
the other hand, the recently observed j22| suppression of high px hadron yields at forward 
rapidities in dAu collisions at RHIC has the behavior anticipated on the basis of CGC 



ideas 
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26|. All these results provide strong evidence for the CGC physics at 



HERA and RHIC. However, more definite conclusions are not possible due to the small 
value of the saturation scale in the kinematical range of HERA and due to the complexity 
present in the description of dAu collisions, where we need to consider the substructure of 
projectile and target as well as the fragmentation of the produced partons. As a direct 
consequence, other models are able to describe the same set of data (See e.g. Refs. 0, 
In order to discriminate between these different models and test the CGC physics, it would 
be very important to consider an alternative search. An ideal place are the future electron- 
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nucleus colliders, which probably could determine whether parton distributions saturate 
and constrain the behavior of the nuclear gluon distribution. This expectation can easily 
be understood if we assume the empirical parameterization Q'^ = x Qq {^)'^, with the 
parameters Ql = 1.0 GeV^, xq = 0.267 x 10^'* and A = 0.253 as in Ref. Q. We can observe 
that, while in the proton case we need very small values of x to obtain large values of Q^, 
in the nuclear case a similar value can be obtained for values of x approximately two orders 
of magnitude greater. Consequently, nuclei are an efficient amplifier of parton densities. 
The parton density that would be accessed in an electron - ion collider would be equivalent 
to that obtained in an electron - proton collider at energies that are at least one order of 
magnitude higher than at HERA 2!^ . 

Recently, an electron - ion collider at RHIC has been proposed in order to explore the 
relevant physics of polarized and unpolarized electron-nucleus collisions |30]. In particular, 
this collider will explore the high density regime of QCD even though its x — Q"^ range will be 
somewhat less extensive than achieved at HERA. For instance, with energies ^/s = 60 — 100 
GeV, one will access x ~ 10^"^ — 10^'^ for ~ 1 — 10 GeV^, respectively. However, the sat- 
uration scale will be approximately 4.0 GeV^ at small x, low and A = 197. Furthermore, 
as pointed in Ref. [3^ , in principle all the inclusive and semi-inclusive observables that were 
studied at HERA can be studied at eRHIC. This collider is expected to have statistics high 
enough to allow for the determination of the logarithmic slopes with respect to x and 
of the total and longitudinal structure functions. In particular, the longitudinal structure 
function is expected to be measured for the first time in the kinematical regime of small x, 
since the electron - ion collider will be able to vary the energies of both the electron and ion 
beams. It will be possible to check predictions made by CGC inspired models (which have 
been extensively tested at HERA) for the behavior of these observables. 

In this paper we study the behavior of the total, longitudinal and charm structure func- 
tions in the kinematical region which will be probed in electron-ion collisions at RHIC 
considering a generalization for nuclear targets of the saturation model proposed by lancu, 
Itakura and Munier (IIM model). Moreover, we estimate the logarithmic slopes of the total 
and longitudinal structure functions at different values of the atomic number. We hope that 



our results contribute to the pL 



Refs. 
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of future eA experiments (for previous studies see 
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We start from the space-time picture of the electron-proton/nuclei processes j48|. In the 
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FIG. 1: The r-dependence of the photon-nucleus overlap functions, normalized by A, for different 
values of the atomic number (x = 10~^ and = 1 GeV^). 

rest frame of the target, the QCD description of DIS at small x can be interpreted as a 
two-step process. The virtual photon (emitted by the incident electron) splits into a qq 
dipole which subsequently interacts with the target. In terms of virtual photon-target cross 
sections ct^.l for the transversely and longitudinally polarized photons, the F2 structure 
function is given by 

F,(x,g^) = ^(a. + a.) (1) 

and 

(Tt,l = j (frdz\^TA'^,z,Q'^)\^ (Jdip^x.r), (2) 

where \l/r,L is the light-cone wave function of the virtual photon and adip is the dipole cross 
section describing the interaction of the qq dipole with the target. In equation Q t is the 
transverse separation of the qq pair and z is the photon momentum fraction carried by the 
quark (For details see e.g. Ref. 0]) . 

The dipole hadron cross section (Jdip contains all information abou^the target and the 
strong interaction physics. In the Color Glass Condensate formalism 
computed in the eikonal approximation and is given by: 



, me laj 



adip can be 



adip{x,r) = 2 J d^bAf{x,r,b) , (3) 
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where is the dipole-target forward scattering amphtude for a given impact parameter b 
which encodes all the information about the hadronic scattering, and thus about the non- 
linear and quantum effects in the hadron wave function. The function JV can be obtained 
by solving the BK (JIMWLK) evolution equation in the rapidity Y = ln(l/2;). It is useful 
to assume that the impact parameter dependence of JV can be factorized as A/'(x, r, b) = 
Af{x, r)S{b), so that <7dip{x, r) = aoAf{x, r), with ctq being a free parameter related to the 



le dipo 



non-perturbative QCD physics. Several models for t 

in the literature in order to fit the HERA data Q, Q, 0, Q, 



e cross section have been used 



19|. Here we will consider 



only the model proposed in Ref. |3| where the dipole-target forward scattering amplitude 
was parametrized as follows, 

1 -exp-''^°'('''^^») , for rQ,(a;) > 2, 

where the expression for rQs{x) > 2 (saturation region) has the correct functional form, as 
obtained either by solving the BK equation [2, 5], or from the theory of the Color Glass 
Condensate [6]. Hereafter, we label the model above by IIM. The coefficients a and b are 
determined from the continuity conditions of the dipole cross section at rQs{x) = 2. The 
coefficients 7^, = 0.63 and k, = 9.9 are fixed from their LO BFKL values and ctq = '^'^Rp, 
where Rp is the proton radius. In our further calculations we shall use the parameters 
Rp = 0.641 fm, A = 0.253, xq = 0.267 x 10"^ and TVq = 0.7, which give the best fit 
result. Rece ntly , this model has also been used in phenomenological studies of vector meson 
production j^| and diffractive pro cesses Isi ] at HERA as well as for the description of the 
longitudinal structure function |52j |. 

Some comments related to IIM model are in order here. Firstly, it is important to 
emphasize that this model is constructed by smoothly interpolating between two limiting 
cases which are analytically under control, but have leading order accuracy. The first line of 
Eq. (0)) is obtained from the solution of the BFKL equation via a saddle point approximation, 
valid for very high energies and very small dipole sizes. The so obtained solution is then 
further expanded under the assumption that the dipole sizes are close to the saturation 
radius. Secondly, it is valid only in a limited range of virtualities, such that the DGLAP 
evolution can be disregarded and the scaling solution of the BFKL equation can be used. 
Moreover, the free parameters in the IIM model have been determined without including 
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FIG. 2: The r-dependence of the photon-nucleus overlap functions, normalized by A, for different 
values of the atomic number (x = 10~^ and = 10 GeV^). 

the charm, i.e. considering only three active flavors. Consequently, our predictions for the 
nuclear charm structure function should be considered as a rough estimate of this observable. 
However, we believe that our main conclusions related to are not modified if a new fit 
including the charm is performed. 

We generalize the JIM model for nuclear collisions assuming the following basic trans- 
formations: ao ^ aQ = X ctq and Ql{x) ^ = A^ x Ql{x). Moreover, in order to 
estimate the contribution of the saturation physics in what follows we present a comparison 
between the full JIM model and the predictions from linear physics, obtained extrapolating 
the expression of A/'(a;, r), valid for rQs{x) < 2, for all kinematical range. In our cal- 
culations the impact parameter dependence of the scattering amplitude, which is mainly 
associated to non-perturbative physics, is disregarded. Basically, following Ref. [l^, we 
shall treat the nucleus as a homogeneous disk of radius Ra- Consequently, with this model 
we cannot discuss the expansion of the transverse size of the target with increasing energy. 
Furthermore, we are assuming that the A ^-dependence of the nuclear saturation scale is 
preserved by the evolution. This assumption is valid in the fixed coupling case, where the 
A scaling of the initial condition survive without change the evolution. On the other hand, 
for a running coupling the A^ scaling holds only in a limited kinematical range, with the 
nuclear dependence becoming weaker at larger energies [l^. In this case, we can expect 
that the A^ dependence of the nuclear saturation scale will be slightly modified at eRHIC. 
It is important to emphasize that more sophisticated generalizations to the nuclear case can 



be used (See e.g. Refs. Consequently, this work should be considered as 

an exploratory study which has as main goal to present a semiquantitative estimate of the 
CGC effects in the future eA collider. In a full calculation we must use the solution of the 
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BK equation, obtained without neglecting the impact parameter dependence as well as an 
initial condition constrained by current experimental lepton-nucleus data. 

Before presenting our results for the total, longitudinal and charm structure function, we 
can investigate the mean dipole size dominating each of these cross sections. We define the 
photon-nucleus overlap function, normalized by the atomic number A, as follows 



where i = T, L characterizes transverse and longitudinal photons. In particular, we also 
calculate the overlap function associated to the total structure function, H2{r,x,Q'^) = 
Ht {r,x,Q'^) + Hl {r,x,Q^), and the overlap function associated to the charm structure 
function, Hc{r,x,Q'^), which is calculated using mj = nic = 1.5 GeV. In Figs. [TJ|21and|Sl 
we show the distinct overlap functions (normalized by A) as a function of the dipole size 
for different values of x, A and Q^. The first aspect that should be emphasized is that 
although the overlap functions have been normalized by A, they are strongly A dependent. 
This behavior is expected when we consider the full prediction of the IIM model. However, 
our results demonstrate that this dependence is also present when we calculate the overlap 
functions using the linear approximation. It is associated to the [Q^ Ai^W"^^ dependence of 
the dipole scattering amplitude, where 7e// = 7s + ^^^^^xy"^ model. For 7e// = 1 

we will have an A 3 -dependence for Af, which combined with the Ai-dependence of ctq implies 
a linear A dependence for the dipole cross section in the linear regime. If normalized by A we 
will obtain an A independent overlap function in the linear regime. However, since 7e// < 1 
in the IIM model, the overlap function (and the corresponding observables) has an A s 
dependence, i.e. it decreases with increasing atomic number. This behavior is observed in 
the figures. In Fig. ^we estimate the overlap functions for fixed x (= 10^^) and (= 1 
GeV^) and two values of the atomic number. In order to illustrate our results we define 
the quantities Vg^A = '^/Qs,a, fs,p = '^/Qs,p and Vq = 2/Q, which are directly associated 
to the nuclear, proton saturation scale and photon virtuality, respectively. The value of 
these quantities is indicated in Figs. ^ and |21 by vertical lines. We have Vg^A < 'fs,p < 
at = 1 GeV^, while vq < rg^A < fs,p at = 10 GeV^. We can see that the charm 
overlap function is peaked at approximately r ^ 0.07 fm, which agrees with the theoretical 



expectation that the cc pair has a typical transverse size ~ 1/yU, where ^ = ^Q"^ + 4m^ 

(See similar discussion in Ref. Therefore, the main contribution to the cross section 
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comes from the small dipole sizes, i. e. from the region where the saturation effects are 
small (linear regime). This expectation is confirmed by the behavior of the overlap function 
H(. which is the same in the linear and full predictions for fixed A. Therefore, we should 
expect that the modifications in the charm structure function due to saturation effects will 
be small. This is expected since the typical scale for charm production, /x^, is larger than 
the saturation scale Ql ^ in all kinematical range of eRHIC, i. e. at this collider it is expected 
that the linear regime dominates the heavy quark production (See discussion in Ref. 21[). 
On the other hand, for light quark production a broader r distribution is obtained, peaked at 
large values of the pair separation, r ^ Vg^p {ts^a) aX A = 1 (197), implying that saturation 
effects contribute significantly in this case. This same feature is observed in the behavior 
of the overlap functions H2 and presented in Fig. ^ In this case we can see that the 
saturation effects suppress the contribution of large dipole size, as expected theoretically. 
Moreover, these effects become more important for smaller values of x and larger A. We 
can observe that the area under the curve is significantly reduced by the saturation effects, 
which implies that the associated observable will be strongly modified by these effects. In 
Fig. |21we present the behavior of H2 and Hl for = 10 GeV^. In this case we observe that 
the distributions peak at smaller values of the dipole size, r ^ Vq, with the contribution 
of large dipole size being reduced. It implies that for A = 1 the full and linear predictions 
are identical. For large nuclei {A = 197) saturation effects still contribute, which implies a 
reduction of the area under the curve and a modification of the associated observable. The 
charm overlap function (not shown) has identical behavior for the linear and full predictions 
and the two values of the atomic number. We can conclude that the saturation effects are 
strongly reduced for large values of Q^. 

In Fig. 121 we present the H2 and Hl overlap functions for two values of x and Q^. In 
this figure we show only the full predictions. As discussed before: (a) increasing the 
distributions peak at small values of the dipole size; and (b) the area under the curve is 
reduced increasing the atomic number. The main aspect of this figure is that it allows 
to analyse the x-dependence of the overlap function. We observe that decreasing x the 
overlap function grows, with the growth being smaller for large nuclei. Consequently, we 
expect that the associated observables increase at small x, with a smaller slope for A = 197. 
These expectations are confirmed in Fig. |3J where we present the x-dependence of the total, 
longitudinal and charm structure functions. We can see that the full and linear predictions 
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FIG. 3: The r-dependence of the photon-nucleus overlap functions, normalized by A, for different 
values of the atomic number, x and Q^. 

for the charm structure function are identical, having the A-dependence characteristic of 
the IIM model. On the other hand, the behavior of and is strongly modified by 
the saturation physics, with the effect decreasing for larger Q^. In order to obtain a more 
precise estimate of the modification in the observables, in Fig. Elwe present the ratios 
and Rpj^ between the full and linear predictions for F2 and Fi, respectively. We consider 
three typical values of the atomic number. As expected, the contribution of the saturation 
physics increases at large nuclei and smaller values of x. In particular, for values of x around 
10^^ we predict a reduction of about 50% in the total and longitudinal structure functions. 
In Fig. ini we present the behavior of the ratio between the nuclear and proton structure 
functions, R{x, Q"^) = F^{x, Q^)/-^^^, Qflj as a function of x and Q^. For comparison, the 



predictions of the EKS parameterization 53|, which is a global fit of the nuclear experimental 
data using the DGLAP equation, is also presented. In particular, in Fig. |B1 (a) we show the 
behavior of this ratio as a function of x at = 2.5 GeV^, assuming two different values of 
A: 197 and 40. We can see that, similarly to EKS parameterization, the ratio obtained using 
the IIM model generalized for nuclear targets (IIMn) decreases when A increases. The main 
difference between the predictions is the behavior of the ratio R at small x. While the EKS 
parameterization predicts that the ratio is constant in this limit, the IIMn one predicts that 
the ratio still decreases at smaller values of x. In Fig. IHl (b) we present the dependence 
of the ratio at x = 10~^. We find that while this behavior in the EKS parameterization is 
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FIG. 4: Nuclear structure functions as a function of x for different values of A and 



directly associated to the DGLAP evolution, in the IIMn prediction it is associated to the 
saturation and geometric scaling regime. We see that the predictions differ significantly at 
small values of Q^, where the saturation physics dominates. 

Other observables of interest to study the CGC physics are the logarithmic slopes of F2 
and Fl with respect to x and Q^. It is mainly motivated by the strict relation between the 
gluon distribution and the scaling violations of the total structure function at leading order 
in the DGLAP formalism [s^. In the dipole formalism the scaling violations are directly 
related to the dipole cross section (See, e.g. Ref. 



X, r 



(6) 



Therefore, this observable can be useful to address the boundary between the linear and 

n 

saturation regimes [3J| . This expectation is easily understood. As this observable is strongly 
dependent on the dipole cross section and it presents distinct behavior for Q"^ > Ql and 
< Q^, its experimental analysis would allow to test the x and A dependence of the 
saturation scale. On the other hand, the logarithmic derivative of F2 with respect to x is 
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directly related with the power of growth of this structure function at small x. Basically, if 
we parameterize the total structure function using F2{x,Q'^) = we obtain that 

i.e. this logarithmic slope is directly related to the effective Pomeron intercept (for a similar 
analysis in ep collisions see, e.g. Ref. [s^). In Fig. Wi^) we present the dependence of the 
effective intercept for x = 10^^ and different values of the atomic number. For comparison, 
the prediction of the GBW model generalized for nuclear targets is also presented. For 
small values of both models predict a similar dependence for A. The main difference 
occurs in the region of large values of where their predictions are not expected to be 
valid. However, an A-dependence for A is observed, A being smaller for large nuclei. This 
behavior can be understood considering the prediction of the IIM model for F.f in the region 
rQs{x) < 2 (linear regime). In this case we obtain that the effective power X{x,Q'^,A) is 
proportional to '^"■^^ ^ lny4/ln(l/x) (plus positive A independent terms), which is negative, 
since 7e// < 1, and grows in modulus with A. Moreover, this agrees with our previous 
results, where we have found that the growth of the nuclear structure function at small x 
decreases at larger A. In El^b) the x dependence of the intercept for = 1 GeV^ is shown. 
In this case we only consider the IIM model and show its linear and full predictions. We can 
see that the linear predictions for A are similar. On the other hand, the saturation effects 
imply that A decreases at small x, its reduction being stronger for A = 197. 

In Fig. ISlwe present the x dependence of the logarithmic derivatives of F2 and Fl with 
respect to Q^. These derivatives, as well as A, have been evaluated numerically using the 
DFRIDR routine j^, which is based on the Richardson's deferred approach to the limit. 
We can see that both derivatives have similar behavior, with the linear and full predictions 
being identical for large x. For the F2 slope and A = 1, the difference between the linear 
and full predictions starts at x ~ 10~^, increasing at smaller values of x. On the other 
hand, at A = 197 both predictions differ at values of x smaller than ^ 10^^, with a large 
difference between the predictions at small x. For the F^ slope and A = 197 we have a 
similar behavior, but with the difference between the predictions starting at x ~ 10~^. At 
y4 = 1 we can see that the full prediction is larger than the linear one for the x range of the 
figure. We have found that at smaller values of x the linear prediction becomes larger than 
the full one. Finally, in Fig. El we present the x dependence of the logarithmic derivatives 
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FIG. 5: Ratio between the full and linear predictions for the different nuclear structure functions. 
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FIG. 6: Ratio R{x, Q^) = F.^{x, Q^)/F^{x, Q^) as a function of (a) x and (b) Q^. The predictions 
of the EKS parameterization are shown for comparison. 



of F2 and Fl with respect to x. We can see that while the linear predictions grow at small 
X, the full predictions present a smaller slope. In particular, at A = 197 these observables 
are almost x-independent when we consider the saturation effects. 

As a summary, in this paper we have studied the predictions of CGC physics for electron 
- ion collisions at high energies, using a generalization for nuclear targets of the lancu- 
Itakura-Munier model which describes the ep HERA quite well. We have estimated the 
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FIG. 7: Effective intercept as a function of: (a) and (b) x. 

nuclear structure function F2{x-,Q'^\ as well as the longitudinal and charm contributions. 
Moreover, we have investigated the behavior of the logarithmic slopes of the total and 
longitudinal structure functions in the kinematical region of the future electron - ion collider 
eRHIC. Our results indicate that the experimental analysis of these observables in the future 
electron - ion collider could discriminate between linear and saturation physics, as well as 
constrain the behavior of the saturation scale. Our analysis was restricted to inclusive 
observables. However, the CGC physics also strongly modifies the behavior of exclusive 
observables, as verified, for instance, in diffractive processes at HERA. Studies of diffractive 
interactions in eA interactions are still scarce. Some examples are those performed in Refs. 
40l |46|, where the diffractive photopro duct ion of heavy quark and vector mesons in eA 
collisions were studied. In a forthcoming publication we calculate the nuclear diffractive 
structure function, finding that the analysis of this observable can be useful to constrain the 
CGC physics. 
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